OS 

o 
o 

(N 



> 
(N 
00 

o 

\o 
o 

a,: 

^— > ■ 
c ■ 

=5 

cr 



X 



Entanglement in general two-mode continuous-variable states: local approach and 

mapping to a two-qubit system 

H.-C. Liifl and A.J. FisheiQ 

UCL Department of Physics and Astronomy and London Centre for Nanotechnology, 
University College London, Cower Street, London WC1E 6BT, U.K. 

We present a new approach to the analysis of entanglement in smooth bipartite continuous- 
variable states. One or both parties perform projective filterings via preliminary measurements 
to determine whether the system is located in some region of space; we study the entanglement 
remaining after filtering. For small regions, a two-mode system can be approximated by a pair 
of qubits and its entanglement fully characterized, even for mixed states. Our approach may be 
extended to any smooth bipartite pure state or two-mode mixed state, leading to natural definitions 
of concurrence and negativity densities. For Gaussian states both these quantities are constant 
throughout configuration space. 



There has been growing interest in the quantification of 
entanglement in quantum systems. However, many sys- 
tems of interest have continuous, rather than discrete, 
degrees of freedom The general characterization of 
entanglement in such continous- variable systems is a very 
difficult problem; most of the known results are for Gaus- 
sian states, where the logarithmic negativity Q can be 
calculated for arbitrary bipartite divisions [3[. The en- 
tanglement of formation [4] is known exactly only for 
symmetric bipartite Gaussian states [j| , and further mea- 
sures of entanglement in multipartite Gaussian states are 
a topic of active current research 1, 6] . For non-Gaussian 
states there is some progress in finding criteria for entan- 
glement 0, but much less in quantifying it. 

In this letter we present a new approach to this prob- 
lem, allowing entanglement to be quantified locally in 
arbitrary (including non-Gaussian) continuous-variable 
states. We concentrate on entanglement localized near 
particular regions in configuration space, which we anal- 
yse via a thought experiment in which the entangled state 
is first measured to localise it. This corresponds to a par- 
ticular type of projective filtering, used to identify the 
distribution of entanglement in a state which has a pre- 
existing bipartite structure. This contrasts with previous 
studies [1, 9|, U3] of the entanglement of a finite region of 
space with the rest of the system. We show that in the 
limit where the sizes of the measured regions become very 
small the description of each mode in the system becomes 
isomorphic to a single qubit. In two-mode systems sim- 
ple expressions result for entanglement monotones (con- 
currence and negativity), yielding natural definitions for 
corresponding densities in configuration space. 

The filtering process. Let Alice and Bob share a state 
of two distinguishable one-dimensional particles (or two 
modes) . Alice can measure only the position of her par- 
ticle or mode (coordinate qa), Bob the position of his 
(coordinate qs)- They filter their state by determin- 
ing whether or not the particles are found in particu- 
lar regions of configuration space, and discard instances 
in which they are not. We refer to the resulting sub- 



ensemble as the 'discarding ensemble'. For example, if 
Alice measures whether her coordinate is in a sub-region 
a, the corresponding projector is 



E„ 



\qA){qA\ dqA ® i 



Bob 



(1) 



where leob is the identity operation for Bob's particle. 
The density matrix po in the discarding ensemble after 
the measurement is 



Pd{o) = 



EapEa 

Pa 



(2) 



where p a = Tv[E a p] is the probability of finding Alice's 
particle in a. Let the entanglement of po be Ed- 

On the other hand if Alice chooses not to discard the 
system when it is not in region a (the 'non-discarding 
ensemble'), the appropriate density operator pnd is 

Pnd = E a pE a + E a ,pE a >, (3) 

where the complementary projector E a i is defined as 



E„> = 1 — E„ 



dq A \qA) (qA\ ® lBob- (4) 



Eq. §3§ describes a mixed state, differing from the orig- 
inal p in that off-diagonal elements connecting q A G a 
and qA ^ a have been set to zero. If all the operators 
available to Alice have support only in region a (i.e. if 
she can neither measure her particle's properties, nor ma- 
nipulate it in any way, except when it is in a) then the 
second component E a ipE a > contains no usable entangle- 
ment. Alice can distinguish the two portions of pnd by 
a local measurement, so the entanglement End in the 
non-discarding ensemble is just End = PaEo- We shall 
focus on calculating Ed in this paper, noting that End 
can be simply obtained from it. 

Pure states: preliminary measurement on Alice's parti- 
cle only. Suppose the state p is pure; so is pr>. It is there- 
fore straightforward to calculate its entanglement from 



2 



the von Neumann entropy of the corresponding reduced 
density matrix p^ D = Tr b[pd]- Suppose further that 
the initial filtering is performed only by Alice, by deter- 
mining whether q A lies in the region q A — a < q A < q A +ci, 
and all instances in which this is not the case are dis- 
carded. Now, since a is to be very small, Alice's original 
reduced density matrix p^ (before the measurement) in 
the neighbourhood of q A can be expanded (provided it is 
smooth in configuration space) as 

P [A) {lA, q' A ) = P (A Hn A +x,q A +y) (5) 
= P (A) oo + P (A) m x + P (A) oiV + P (A \iXy + 0(x 2 ,y 2 ), 



where 



P [A) 

r nm 



^wr p[A){qA ' q ' A K= q ' A =- qA - ^ 



Within region a, p { > D is obtained by rescaling p { ' ac- 
cording to ©, where Ti[E a p] = 2a[p ( - A \ + 0(a 2 )]. 

Now seek right eigenfunctions 4> n of p( A ) D within the 
region a. Expanding <f> n as a power series 



4> n (x) = a n + b n x + 0(x 2 ), 



(7) 



the eigenfunction condition becomes a matrix-vector 
equation operating on the expansion coefficients 
{a, b, . . .}. To order a 2 , the non-zero eigenvalues are: 



Ai 
A 2 



1-Ao 



3^oo (P 



(A) JA) 
llP 00 



P oiP ioJ- 



(8) 



So to lowest order (a 2 ) the eigenvalues, and hence the 
von Neumann entropy, of p^ D are entirely determined 
by the quantity e = A 2 ; the corrections due to higher 
eigenvalues, or to higher-order terms in Eq. ((5]), affect 
the result only to order a 4 . Specifically, the von Neumann 
entropy is 



S = h(e) = -[elog 2 (e) + (1 - e)log 2 (l - e) 



(9) 



Note that if Alice's state is pure, p^ A \iP^ A \a = 
P oi P id so & i s zer0 as we would expect. 

Pure states: preliminary measurement on both par- 
ticles. It is possible to generalise this analysis to the 
case where both Alice and Bob make a prelminary mea- 
surement to localise their particles, within regions {q A '■ 
q A -a <q A < q A +a} and {q B :q B -b<q B < q B + b} re- 
spectively In that case one can expand p as a joint power 
series in {q Al q' A , q B , q' B }, calculate the reduced density 
matrix p( A ' (also as a power-series expansion) and pro- 
ceed as above. However, further insight can be obtained 
by an alternative approach. Define for both Alice and 
Bob two-dimensional state spaces consisting of 



4>Ao(xi) - 
4>Bo(x 2 ) 




&Al(xi) z 
4>Bl(x 2 ) 




which are orthonormal on the intervals —a < x\ < a 
and —6 < x 2 < b respectively; (f>o represents the constant 
component of the wave function, and 4>i the spatially 
varying part. So long as terms varying as x 2 or higher 
can be neglected, a Taylor expansion of the joint state to 
linear order ([5]) is equivalent to expanding tp in the basis 
(flC)]) , thereby reducing the joint system to a two-qubit 
one. It can be shown that the third largest eigenvalue of 



P 



(A) 



(corresponding to the extent to which the two-qubit 
approximation fails) is now of order (ab) 4 . 

We can now use any of the standard measures of the 
entanglement of the two-qubit system. For pure states, 
the tangle [ll| is 



1 



ab\ipoo\ 2 



Aa 2 b 



2u2 



■(^ooV'ii - V>oiW 



where 

Ipr, 



3q A n dq B m ^ A=1 iA,qB=q B - 



(11) 



(12) 



The prefactor in (jllj) comes from the normalization con- 
dition J^ a dxi J. dx 2 \fp(xi, x 2 )\ 2 = 1. The entangle- 
ment is therefore h ((1 - Vl~7)/2) = h (t/4 + 0(r 2 )) . 
By analogy with the definition of concurrence C = y/r 
for two-qubit states [HJ , we define a concurrence density 
c = C / (ab) such that r — (cab) 2 ; then 



3/9, 



■[piiooPoon + PooooPnn 



0000 



"PioooPom — PoiooPion 



1/2 



(13) 



where 

Pijkl 



Qi QJ Qk Ql 

^^^^ p{qA ^ qB;qA ^ 



(14) 



for pure states. 

The negativity M of the filtered state can also be com- 
puted; we define this as the sum of the magnitudes of 
the negative eigenvalues A^ of the partially transposed 
density matrix p Tu , 



E 



I A, 



(15) 



iS.t. Ai<0 



For pure states negatitivity and concurrence are simply 
related (HQ: N = C/2. 

The accuracy of the two-qubit approximation is guar- 
anteed (for sufficiently small a and b) by the fact that 
each party's reduced density matrix has only two non- 
zero eigenvalues of the density matrix to order (ab) 2 . 

Mixed states. The mapping to a two-qubit system ap- 
plies also to mixed states, where exact recipes for the 
entanglement of formation [l2| and other measures are 
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known. We find all eigenvalues fa of pp (as defined in 
[lH) are at leading order proportional to (ab) 2 so again 
there is a well-defined concurrence density c. These lead- 
ing terms (and hence the concurrence C) can be found 
by solving a quartic, although its roots are not simple in 
general. 

Particularly simple expressions can be found for the 
local negativity of a mixed state. We define the quantities 



D 2 

Ci 



(pnooPoooo — PoiooPiooo)/(3poooo); 

(PoonPoooo — PoooiPooio)/(3poooo); 

[(pooooPoioi — PoooiPoioo)pioio (16) 
+Pooii(poiooPiooo — P0000P1100) 
+Pooio(poooiPiioo — poioipiooo)] /(9poooo)- 



If C\ > 0, the negativity depends on a and b individually, 
not just on the product ab. If we seek the maximum neg- 
ativity while keeping ab fixed, we find this occurs when 
a 2 D\ — b 2 D 2 , and obtain a corresponding negativity den- 
sity n = 7V max /(a6) given by 



y/Ci + D 1 D 2 - ^D X D 2 



(17) 



For pure states, D\ and D 2 vanish, and we recover the 
pure-state negativity expected from equation (fT3|) . If 
C\ < there is no negativity to order (ab) ; the full treat- 
ment of the negativity will be discussed in a subsequent 
paper. 

Example: Entanglement of Gaussian states. The char- 
acteristic function C is defined in terms of the Weyl op- 
erator C (taking fi = 1) through 

C(X,P) = Tr(pC(X,P)); C(X,P) = e i{x « +p§ \ 

(18) 

where the position operator is denoted by q and the mo- 
mentum operator by p. A state p is said to be Gaus- 
sian when its characteristic function is a Gaussian in 
phase space. This important set of states includes both 
thermal and 'squeezed' states of harmonic systems and 
plays a key role in several fields of theoretical and ex- 
perimental physics; we use them as an example because 
their entanglement properties are better understood than 
those of other continuous-variable systems, while recog- 
nising that our approach is general. The corresponding 
configuration-space density matrix can be written (l5j 



P{q\l') = Ciexp[-g T Lg - q' T Lq' 



(19) 



--( q -q>fM(q-q>) 



(q-q') T K(q + q') 



where Ci is a normalization constant, and where L, M 
and K are real N x N matrices (for N modes) with L 
and M symmetric, while K is arbitrary. For thermal 
states of two similar but distinguishable oscillators each 
having mass m, angular frequency uj and coupling spring 
constant K with corresponding dimensionless coupling 




« \ q=0 



q=2.4 % 



(a) 



(b) 



FIG. 1: Entanglement properties as a function of region size 
for a Gaussian ground state with a = 10 and m = uj — 1, 
in the case where both Alice and Bob make preliminary mea- 
surements and the reigion sizes are chosen to be the same, (a) 
Entanglement S (dimensionless) as a function of region size 
2a (in units of (mo;) -1 ^ 2 ; log scale) for two different positions 
(data points); the entanglements contained in the effective 
two-level systems constructed from the two largest eigenval- 
ues of p^ A ^ are also shown (dashed lines), (b) The concurrence 
density c (in units of mui), computed from the entanglement 
of formation by inverting the relation S = h((cab) 2 /4); note 
how it saturates to the exact result predicted by (I20|> (hori- 
zontal line) for small regions. 



a = 2K jmiJK K = and the values of L and M arc 
given by 0,[lB|. We adopt (mw)" 1 ' 2 as our length unit. 
First we consider the ground state. Eq. ([TB"]) becomes 



V2r 



G 



1 + 2a - Vl + 4a. 



(20) 



So, for a given small a and &, the entanglement depends 
only on a and the fundamental length unit; it is indepen- 
dent of the location of the centres (JIa^b) of the mea- 
surement regions. The lack of dependence on {JiajQb) 
is a special feature of Gaussian states, and this result 
does not hold in general; note also that the concurrence 
density can be made arbitrarily large by increasing a. 

Figure Q] shows the variation of the entanglement S 
(calculated numerically) with the region size. Note how 
the entanglement saturates to the full value given in [lij 
for large regions, while for small regions it reduces to the 
value predicted by ([217)1 . To obtain entanglements of a 
substantial fraction of one ebit, it is necessary to choose 
a region size comparable to the fundamental length unit 
of the oscillator; around this point the two-qubit approx- 
imation is just starting to break down. Calculational 
details and further results are described in [17|. 

For mixed Gaussian states, we find that the concur- 
rence density and negativity density are again indepen- 
dent of position. In Figure [2] we plot both quantities as 
a function of temperature for thermal states of two os- 
cillators with different coupling strengths. We also show 
the conventional (global) negativity for this state; note 
that both local entanglement densities vanish at the same 
temperature as the global negativity, showing that for 
this set of states, those which are entangled from the 
global point of view are also entangled by our local mea- 
sures. It would be very interesting to know whether 
this property applies more generally. We also note that 
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FIG. 2: Negativity density n (full lines; in units of mui ) 
and concurrence density c (long dashes; in units of mui) as 
a function of temperature T (in units of u>) for a thermal 
state of the two-oscillator system discussed in the text having 
kg = m = u) = 1, and for two different values of the coupling 
a. The global negativity M g [|| is shown for comparison (short 
dashed lines; dimensionless) . 



the concurrence density is exactly twice the negativity 
density; this is a special property of two-mode Gaussian 
states, which will be discussed further elsewhere. 

Extraction of the local entanglement. Methods of ex- 
tracting the entanglement from a squeezed continuous- 
variable state into a p air of two-level system were pre- 
viously studied in [3 S HI Using our mapping to an 
effective two-qubit system, we can swap the entangle- 
ment in a small region of the continuous wavefunction to 
local qubits (i.e. true two- level systems). Remembering 
that the states 4>o an d (pi drop to zero outside the re- 
gion [—a, a] we find that the Pauli operators X and Y 
of the effective qubit can be represented in terms of the 
canonical position and momentum operators q and p by 



it 

X = q; 

a 



Y 



2a 



p- 



(21) 



The experiment could be performed as follows: first, lo- 
calize the continuous degree of freedom (for example, 
through a homodyne measurement in the case of an elec- 
tromagnetic field mode), then perform a SWAP opera- 
tion by composing three controlled-X gates: 



U s 



WAP 



exp 



j(*„-U(^«-l) 



x exp 



x exp 



■J(0"x-1)(-^P-1) 



ITT 



4 ( ( 7,-l)( — q-1) 



,V3, 
a 



(22) 



where {<j x ,<j v ) are Pauli operators for a local qubit. Us- 
ing this procedure one could therefore extract the full 
two-level entanglement shown in the dashed curves of 
Figure [2a ) provided the swap operation is successful. 

Discussion. After filtering in configuration space as 
we have described, it is simple to characterize the en- 
tanglement in a continuous-variable system. One could 



generalize our results to multimode systems by making 
multivariate Taylor expansions of p or by using the iso- 
morphism of each mode to one qubit. Our approach 
leads naturally to quantities characterizing the entangle- 
ment that scale with the extent of the measurement re- 
gion used in the filtering, and hence to the identification 
of well-defined densities. However, the concurrence and 
negativity are not extensive, in the sense that the sum of 
these quantities over all the sub-regions of configuration 
space does not yield the full entanglement of the original 
system. This point is discussed further in [171]. It will 
be interesting to characterize further the relationship be- 
tween the local and global views of continuous-variable 
entanglement; in any case our results open a wide range 
of non-Gaussian states to further study. 
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